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Polarization and radiation reaction (RR) effects in the interaction of a superintense laser pulse (/ > 
10 23 W cm ) with a thin plasma foil are investigated with three dimensional particle-in-cell (PIC) 
simulations. For a linearly polarized laser pulse, strong anisotropies such as the formation of two 
high-energy clumps in the plane perpendicular to the propagation direction and significant radiation 
reactions effects are observed. On the contrary, neither anisotropies nor significant radiation reaction 
effects are observed using circularly polarized laser pulses, for which the maximum ion energy exceeds 
the value obtained in simulations of lower dimensionality. The dynamical bending of the initially 
flat plasma foil leads to the self-formation of a quasi-parabolic shell that focuses the impinging laser 
pulse strongly increasing its energy and momentum densities. 

PACS numbers: 52.38.Kd, 52.65.Rr 



The radiation pressure generated by ultraintense laser 
pulses may drive strong acceleration of dense matter, as 
experimentally shown in various regimes Radiation 
pressure may thus be an effective mechanism for the gen- 
eration of high-energy ions, especially in the regime of 
extremely high intensities and relativistic ion energies as 
foreseen with the ELI project. In the case of solid-density 
thin foil targets, pioneering particle-in-cell (PIC) simula- 
tions showed that at intensities exceeding 10 23 Wcm~ 2 
and for linear polarization of the laser pulse radiation 
pressure dominates the acceleration yielding linear scal- 
ing with the laser pulse intensity, high efficiency and 
quasi-monoenergetic features in the ion energy spectrum 
0| . More recent two-dimensional (2D) simulations for a 
small disk target suggested a potentially "unlimited" en- 
ergy gain for the fraction of ions that get phase-locked 
with the laser pulse 0] ■ 

The above mentioned studies showed that the radia- 
tion pressure dominant acceleration (RPDA) regime is 
very appealing as a route to the generation of relativistic 
ions, but leave several theoretical issues open. Firstly, 
transverse instabilities Q and multi-dimensional effects 
may play a crucial role as shown by 2D simulations Q. 
Secondly, the use of circular polarization (CP) instead of 
linear polarization (LP) quenches the generation of high- 
energy electrons Q , allowing radiation pressure to domi- 
nate even at intensities below 10 23 W cm -2 and leading to 
efficient acceleration of ultrathin foils @ ; it has not been 
shown yet whether the use of CP is advantageous also at 
ultra-high intensities (/ > 10 23 Wcm -2 ), i.e. when the 
radiation pressure generated by the laser pulse becomes 
the dominant mechanism of acceleration both for CP and 
LP. Finally, it has been shown by one-dimensional (ID) 
simulations that radiation reaction (RR) effects may sig- 
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nificantly affect the dynamics of radiation pressure ac- 
celeration both for thick Q and thin targets H, Q , and 
also depend strongly on the polarization [8| . All of these 
phenomena may be affected by the dimensionality of the 
problem, and a fully three-dimensional (3D) approach is 
ultimately needed because e.g. in 2D simulations and for 
LP the laser-plasma coupling is different for S- and P- 
polarization (i.e. for the electric field of the laser pulse 
either perpendicular or parallel to the simulation plane, 
respectively) and the constraint of the conservation of an- 
gular momentum carried by CP pulses holds in 3D only. 

In this paper, we address the role of polarization and 
RR effects in the RPDA regime using fully 3D PIC sim- 
ulations. To our knowledge, these are the first 3D simu- 
lations of ion acceleration in the RPDA regime with RR 
effects included and among the largest and most accurate 
3D simulations reported so far. Our results show that 
even in the RPDA regime CP leads to higher ion energies 
and better collimation than LP, for which an anisotropic 
ion distribution is observed. It is also found that the 
bending of the foil leads to a self-generated parabolic shell 
that focuses the impinging pulse down to an almost A 3 
volume and that the energy density at the focus largely 
exceeds the initial peak energy density. Compared to 2D 
simulations with analogous parameters, the pulse focus- 
ing effect is remarkably enhanced and the cut-off energy 
of ions is increased. RR effects on the ion spectrum are 
found to be negligible for CP but quite relevant for LP 
where they increase the energy cut-off. 

Our approach is based on the numerical solution of ki- 
netic equations for the phase space distributions of elec- 
trons and ions, where RR is included in the motion of 
electrons via the Landau-Lifshitz (LL) force [To[ • Details 
of our RR modeling and numerical implementation in a 
PIC code are given in Refs.{§,[ll|. The effective equation 
of motion for electrons, after neglecting terms which are 
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FIG. 1. Spatial distributions of ions (upper row) and electrons (lower row) at t = 20 T and in the region (\y\, \z\) < 5.7A, for LP 
without (a) and with (b) RR and for CP without (c) and with (d) RR. Ions and electrons are divided into seven populations 
according to their kinetic energy, with the color-bar reporting the lower bound of the energy interval. In the LP case (frames 
(a,b)), the polarization is along the y axis. 
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where r c = e 2 /m e c 2 . The RR force contribution is im- 
portant for ultrarelativistic electrons and is usually dom- 
inated by the first term, while the second term ensures 
the on-shell condition @ . Notice that the dominant term 
has almost the same form also in different approaches to 
RR modeling @, . For a plane wave propagating along 
the x-axis, the RR force is maximum or zero for counter- 
propagating (v x /c — > —1) or copropagating (y x /c — > +1) 
electrons, respectively. 

In order to clarify the new qualitative features due to 
RR effects, we recall that the phase space volume ele- 
ment J evolves according to dJ/dt = JV P -in- It has 
been shown that V p ■ ffi < and therefore the RR 
force leads to a contraction of the available phase space 
volume. The physical interpretation of this property is 
that the RR force acts as a cooling mechanism for the sys- 
tem accounting for the emission of high-energy photons. 
These photons are assumed to escape from the plasma 
freely, carrying away energy and entropy (Tlj . 

We present a total of four 3D simulations each with 
the same physical and numerical parameters but differ- 
ent polarization, with and without RR effects. In these 
simulations, the laser field amplitude has a sin 2 -function 



longitudinal profile with 9A FWHM (where A = 0.8 fj,m 
is the laser wavelength) while the transverse radial pro- 
file is Gaussian with 10A FWHM and the laser pulse 
front reaches the edge of the plasma foil at t = 0. The 
peak intensity at the focus is I = 1.7 x 10 23 Wcm~ 2 
which corresponds to a normalized amplitude oq = 280 
for LP and ao = 198 for CP. The target is a plasma 
foil of electrons and protons with uniform initial density 
no = 64n c (where n c = 7rm e c 2 /e 2 A 2 is the critical den- 
sity), thickness £ = 1A and initially located in the region 
10A < x < 11A. The density n ~ 1.1 x 10 23 enr 3 is 
slightly lower than that of solid targets but the areal den- 
sity n§l has fully realistic values. Moreover, laser pulses 
of ultrahigh contrast arc now available [T^, and refer- 
ences therein] to avoid early plasma formation effects by 
the prepulse, thus a thin plasma with step-like profile is 
not an unrealistic assumption. 

The simulation grid is 1320 x 896 x 896 and the spatial 
step is A/44 for each direction. The timestep is T/100 
where T = X/c = 2.67 fs is the laser period. We use 216 
particles per cell for each species and the total number 
of particles is 1.526 x 10 10 . The runs were performed 
using 1024 processors each one equipped with 1.7 GB of 
memory of the IBM-SP6 cluster at the CINECA super- 
computing facility in Bologna, Italy. 

FigfJ] shows the ion and the electron 3D spatial distri- 
butions at t — 20 T for the LP case without (a) and with 
(b) RR and for the CP case without (c) and with (d) RR. 
The color corresponds to the range in kinetic energy. For 
CP, the ion spatial distribution follows the spatial inten- 
sity profile of the initial laser pulse, has rotational sym- 
metry around the central axis, and a distribution in en- 
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FIG. 2. (x,y) and (x,z) sections of the 3D simulations of the laser pulse-foil interaction [al-3), bl-3), cl-3), dl-3)] and 2D 
simulations for CP and without RR with the same parameters as the 3D simulations [e), f)], all at t = 20 T. Each frame reports 
the color contours of \/E 2 + B 2 (normalized units) in the x — y plane at z — (al-3 and cl-3), in the x — z plane at y = (bl-3 
and dl-3) and in the simulation plane e), f) in the 2D case. Line contours of the ion and electron densities are superimposed 
in the upper and lower frames, respectively. The CP case with RR is almost identical to the CP case without RR and it is not 
reported. 



crgy monotonically decreasing with increasing radial dis- 
tance. The most energetic ions are located near the axis. 
The number of ions having energy £ > 1100 MeV and 
£ > 800 MeV are 2.3 x 10 10 and 9.4 x 10 10 , respectively. 
The electron spatial distribution has an helicoidal shape 
with step A Fig[T] c),d). RR effects play a minor role 
for CP affecting only a small fraction of ultra-rclativistic 
electrons (mostly removing fast electrons behind the foil 
with almost no influence on the ion distribution as seen 
by comparing frames c) and d) in FigO}. 

The (x, y) and (x, z) sections of the total electromag- 
netic energy density and of the ion and electron densities 
for the CP case in Fig|2]a3)-b3) and c3)-d3) evince a self- 
generated parabolic shell wrapping the laser pulse and 
focusing it up to nearly a A 3 volume, so that both the 
energy and the momentum densities at the focus reach 
values more than eight times their peak value in the ini- 
tial laser pulse. This effect is much weaker in 2D simula- 
tions with the same parameters as shown in Figs|5]e)-f). 
Along the axis, the peak value and width of the ion den- 
sity profile are ~ 10?i c and ~ 0.5 A, showing a strong 
rarefaction due to the transverse expansion, potentially 
leading to enhanced acceleration as described in Ref.0|. 

For LP, the peak ion energy is lower than for CP, the 
ion distribution is anisotropic and RR effects are much 
stronger. The most energetic ions (800-1100 MeV) arc 
grouped into two off-axis clumps lengthened and aligned 
along the polarization direction, and their number is in- 
creased in the case with RR as seen by the comparison 
of FigsUJa) and b) and also in Figs[2]al),a2) and bl),b2) 
where sections of the ion density in the (x,y) and (x, z) 
planes are shown. The contours of the electromagnetic 
(EM) energy density in Figs[2] al),a2) and bl),b2) show 
that near the axis most of the laser pulse has been trans- 
mitted through the target. The increased bunching and 
higher density observed in the case with RR may be re- 



lated to the higher ion energies since the local increase 
of the density and therefore of the reflectivity leads to a 
longer and more efficient RPDA phase. This is consistent 
with observing in Figs|3]bl) and b2) that the EM energy 
density is higher behind the two high-density clumps, 
which correspond to the most energetic ions and are sim- 
ilar to the ion lobes observed in Ref.[l]| at lower inten- 
sity and in a regime of strong pulse penetration through 
the foil. The pulse focusing effect by the self-generated 
parabola is present also in the LP case although weaker 
than in the CP case, and presumably reduced as the laser 
pulse breaks through the parabolic shell. 

The differences between CP and LP in the accelera- 
tion dynamics are well explained, for planar geometry 
and non-relativistic ions, by the absence of the oscil- 
lating component of the J x B force for CP which 
maximizes the effect of the radiation pressure strongly 
suppressing the electron heating. The absence of the os- 
cillating component of the J x B force also accounts for 
the very different RR effects. For CP, a steady push of 
the foil and weak pulse penetration are observed, most 
of the electrons move coherently with the foil and in the 
same direction as the laser pulse so that the RR force be- 
comes very small in accordance with Eq. (p} since the elec- 
trons effectively copropagatc with the laser pulse (see also 
HQ). For LP, the J x B-drivcn oscillations allow elec- 
trons to collide with the countcrpropagating laser pulse 
twice per cycle producing temporal maxima in the RR 
force in agreement with Eq.([T]). Our present results indi- 
cate that CP leads to more efficient acceleration, produc- 
ing higher energy and collimated ion beams, and making 
RR effects negligible also in the 3D case, accounting for 
target bending and pulse focusing effects, and for rela- 
tivistic ions. 

Large-scale 3D PIC simulations are limited by the size 
and availability of computational resources both in the 
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number of runs that may be performed and in the achiev- 
able numerical resolution. This last issue may raise 
doubts on the accuracy of 3D results. To gain confi- 
dence on this side, as well as to compare the 3D results 
to those obtained in lower dimensionality, we performed 
2D simulations both with numerical parameters similar 
to those of 3D runs and with higher resolution. The ef- 
fect of increasing resolution and particle number on the 
ion spectra in 2D simulations is shown in Figj3] where 
2D results are reported for the three different polariza- 
tion cases (CP, LP-S 1 and LP-P) and compared with the 
3D results for both LP and CP. The spectra are normal- 
ized to unity both in the 2D and 3D cases. In the CP 
case both numerical and RR effects on the spectrum are 
smaller while in the P-polarization case these effects are 
larger. Changing the spatial resolution from A/44 to A/80 
and increasing the number of particles-per-cell (ppc) for 
each species from 256 to 625 shifts the energy cut-off by 
- 2% in the CP case and by - 15% (~ 20%) in the P- 
polarization case without RR (with RR). The stronger 
effect of the inclusion of RR for the higher resolution 
case may be explained by noticing that RR mostly affects 
the highest energy electrons [8| , which are located in the 
high-energy tail of the distribution function that needs 
a very large number of particles to be resolved prop- 
erly. Nevertheless, the limited resolution does not quali- 
tatively affect prominent features in ion spectra, such as 
the higher ion energy for CP and the relevance of RR ef- 
fects for LP only, leading for this latter case to an higher 
energy of ions with respect to the case of no RR as ob- 
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FIG. 3. Ion spectra from 2D [a)-c)] and 3D [d)-e)] simulations 
with same physical parameters, all at t = 20T. The 2D spec- 
tra are reported for circular (CP, frame a)) and linear (LP) 
"5" (frame b)) and "P" (frame c)) polarization cases. In each 
plot, the blue and red curves correspond to simulations with- 
out and with radiation reaction (RR) effects, respectively. In 
the upper plots of frames a)-c) the numerical resolution (num- 
ber of particles per cell and of points per wavelength) is simi- 
lar to those of the 3D simulations in d)-e), while in the lower 
plots the results for higher resolution are shown. 



served in ID simulations [8j. As a novel feature of 2D 
simulations, P-polarization leads to much stronger RR 
effects than S'-polarization. In fact, for P-polarization 
the electric field can drag a large fraction of electrons 
out in vacuum and towards the laser pulse as the plasma 
foil begins to bend, enhancing the RR effect. 

For CP, the maximum ion energy is higher in the 3D 
case (~ 1600 MeV) than in the 2D case (~ 1400 MeV). 
In turn, the latter value is higher than what found in ID, 
plane-wave simulations, for which we find a broad spec- 
tral peak which at t = 20 extends up to ~ 1100 MeV 
and is centered around a value of ~ 870 MeV. The lat- 
ter value corresponds to the energy Els = {l — l)w p c 2 , 
where 7 = 1 A/1 — P 2 and f3 are obtained from the "light 
sail" model [14j, and references therein] by numerically in- 
tegrating the ID equations of motion for the foil 



d-y/3 _ 2I(t-X/c) (1-/3 
dt nolrripC 2 \ 1 + p 



dX 
~dt 



= 0c. (2) 



From the ID modeling we also evaluate a final ion energy 
for the spectral peak of ~ 1700 MeV that is reached at 
t ~ 90. The 3D simulations could not be extended up to 
the end of the acceleration stage (estimated up to ~ 90 T 
in ID simulations) but, since the efficiency of RPDA in- 
creases with the foil velocity, the energy gain is expected 
to be even larger at the later times (provided that 3D 
effects do not cause an early stop of the acceleration). 
Hence, the comparison at t = 20 shows an overall in- 
crease of the ion energy in the 3D case with respect to 
ID and 2D cases. Part of the energy enhancement can be 
attributed to the reduction of the areal density UqI due 
to the transverse expansion, as was noticed in 2D simu- 
lations supporting the model of "unlimited" acceleration 
. An additional contribution may come from the above 
described focusing of the laser pulse by the deformed foil, 
which is stronger in 3D geometry. This latter effect was 
absent in the simulations of Ref. |3| because a target with 
a radius smaller than the pulse waist was considered. 

In conclusion, with three-dimensional particle-in-cell 
simulations of ultra intense laser interaction with solid- 
density foils we showed that circular polarization im- 
proves ion acceleration also in the radiation pressure 
dominant regime, confirming and extending previous re- 
sults obtained for lower intensity and/or lower dimension- 
ality. In detail, circular polarization leads to the highest 
ion energies, to a symmetrical and collimated distribu- 
tion, and to negligible effects of radiation reaction. In 
addition, the maximum energy of ions in 3D is larger than 
observed in corresponding ID and 2D simulations. This 
enhancement is attributed both to the density decrease 
in the target, as noticed in the "unlimited acceleration" 
model 0] , and to the strong focusing of the laser pulse by 
the parabolically deformed foil. In the linear polarization 
case, lower maximum energies are achieved, the most en- 
ergetic ions are grouped anisotropically into two off-axis 
clumps and radiation reaction effects significantly affect 
the energy spectrum. We expect these findings to be of 
relevance for the design of future experiments on laser 
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